Exact Analysis of Soliton Dynamics in Spinor Bose-Einstein Condensates 
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We propose an integrable model of a multicomponent spinor Bose-Einstein condensate in one 
dimension, which allows an exact description of the dynamics of bright solitons with spin degrees of 
freedom. We consider specifically an atomic condensate in the F = 1 hyperfine state confined by an 
optical dipole trap. When the mean-field interaction is attractive (co < 0) and the spin-exchange 
interaction of a spinor condensate is ferromagnetic (C2 < 0), we prove that the system possesses a 
completely integrable point leading to the existence of multiple bright solitons. By applying results 
from the inverse scattering method, we analyze a collision law for two-soliton solutions and find that 
the dynamics can be explained in terms of the spin precession. 
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In 2002, two groups 0, reported matter- wave soli- 
tons of an atomic Bose-Einstein condensate (BEC). They 
prepared BECs of 7 Li atoms in a small region of an elon- 
gated optical dipole trap, which is an analog of a waveg- 
uide for microwaves. After tuning the strength of inter- 
action between the atoms to a sufficiently large negative 
value, they set the condensate free along the waveguide. 
The solitary wave packets were formed and propagated 
in the guide nondispersively. It is well known that the 
Gross-Pitaevskii (GP) equation with attractive interac- 
tions in a one-dimensional (ID) space, which is also called 
the self-focusing nonlinear Schrodinger (NLS) equation, 
has bright soliton solutions Therefore, they con- 

cluded that the dynamics of the system is actually ID 
so that the matter-wave solitons can be observed. 

Matter-wave solitons in a new field of atom optics Q 
are expected to be useful for applications in atom 
laser, atom interferometry, and coherent atom transport. 
Moreover, it could contribute to the realization of quan- 
tum information processing or computation. When one 
explores these future applications, atomic BECs have an- 
other advantage. That is, atoms have many internal de- 
grees of freedom liberated under an optical trap giv- 
ing rise to a multiplicity of signals. The properties of 
BECs with spin degrees of freedom were investigated by 
many researchers 0, IE • 

In this Letter, we combine these two fascinating prop- 
erties, matter- wave soliton and internal degrees of free- 
dom. We consider BECs of alkali atoms in the F = 1 
hyperfine state, such as 7 Li, 87 Rb, and 23 Na, confined in 
the ID space by purely optical means. Under no external 
magnetic fields, their three internal states nip = 1, 0, — 1, 
where nip is the magnetic quantum number, are degener- 
ate. The dynamics of the spinor condensates is described 
by the multicomponent GP equations within the mean- 
field approximation. Those coupled equations have non- 
trivial nonlinear terms reflecting the SU(2) symmetry of 
the spins. 



When the mean-field interaction is attractive and the 
spin-exchange interaction is ferromagnetic, we show that 
this system possesses a completely integrable point. By 
considering a reduction of results from the inverse scat- 
tering method, we present for the first time the exact 
multiple bright soliton solutions for the system with the 
spin-exchange interaction. The spin-exchange interac- 
tion, which is absent in the systems of Ref. 0, Q be- 
cause of frozen spin degrees of freedom under additional 
magnetic fields, gives rise to spin mixing within conden- 
sates || during soliton-soliton collisions. We analyze the 
collision law for two-soliton solutions and find that the 
soliton dynamics can be explained in terms of the spin 
precession. 

The assembly of atoms in the F = 1 state is char- 
acterized by a vectorial order parameter: <f>(x,t) = 
t), &o(x, t), $_i(x, t)] T with the components sub- 
ject to the hyperfine spin space. The normalization is 
imposed as J dx 4>(x, t)> ■ <f>(x, t) = Nt, where Np is the 
total number of atoms. Here we assume that the system 
is one dimensional: the trap is elongated in the x direc- 
tion such that the transverse spatial degrees of freedom 
are factorized from the longitudinal and all the hyperfine 
states are in the transverse ground state. This quasi- 
one-dimensional regime is achievable . The interaction 
between atoms in the F = 1 hyperfine state is given by 
V{x\ — x 2 ) = 5(xi — X2)(co + C2Fi- F 2 ) , where Fj are the 
angular momentum of two atoms . In this expression, 
Co = (So + 2g;»)/3, c.2 — (#2 — <?o)/3, with the effective ID 
couplings, jlCj 
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where aj are the s-wave scattering lengths in the total 
hyperfine spin / channel, a±_ is the size of the transverse 
ground states, m is the atomic mass, and C = —((1/2) = 
1.4603 • • •. Then, the Gross-Pitaevskii energy functional 
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is expressed as 



h 
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where repeated subscripts {a, (3, a', f3' — 1, 0,-1} should 
be summed up and / = (f x , f v , f z ) T with /' being 3x3 
spin-1 matrices. 

The time evolution of spinor condensate wave func- 
tion t) can be derived from the variational principle: 
ihd t ® a (z,t) = 6E GP /5$* a {x,t). Substituting Eq. © into 
this, we obtain a set of equations: 



ifi&Sii = -^^$ ±1 + (c + c 2 )(|$ ±1 | 2 + |$ | 2 )$ ±1 
2m 

+ (CO -C 2 )|$ T l| 2 $ ± l +P 2 ^ 1 $o> 
?i 2 

i?i9 t $o = -t— 9 2 $o + c |4'o| 2 $o + (co + C2)(|$i| 2 

Z77J 



|$_i| 2 )$ + 2c 2 $5$i*-i. 



(3) 



In this Letter, we consider the system with the coupling 
constants cq — c 2 = — c < 0, equivalently 2go — — g 2 > 0. 
The effective interactions between atoms in a BEC have 
been tuned with a Feshbach resonance In spinor 

BECs, however, we should extend this to alternative 
techniques such as an optically induced Feshbach reso- 
nance |l2| or a confinement induced resonance |lCj , which 
do not affect the rotational symmetry of the internal spin 
states. In the latter, the above condition is surely ob- 
tained by setting a± = 3Caoa 2 /(2ao + a 2 ) in Eq. (1) 
when ao > a 2 > or a 2 > > ao- Recently, such strong 
transversal confinement has been realized in a 2D optical 
lattice where a± ~ tens of nm . 

Introducing the dimensionless form, <fr — > (</>i,\/20o, 
</>_i) T , where time and length are measured in units of 
t — ha±/c and x = h\/a±j2rnc, respectively, we can 
rewrite Eqs. 10 as a 2 x 2 matrix version of the NLS 
equation: 



id t Q + &iQ + 2QQtQ = 0, Q = 
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Since the matrix NLS equation (0J is integrable , the 
dynamical problems of this system can be solved exactly. 
The embedding (J3J and its first application to an atomic 
system with the spin-exchange interaction are the main 
idea of this Letter. 

We remark that a different choice of Q [lj] gives rise to 
coupled NLS equations known as the Manakov model 0] 
which is widely used to describe the interaction among 
the modes in nonlinear optics. 

The general 7V-soliton solution of Eq. JIJ is obtained 
through a reduction of a formula derived by the inverse 



scattering method (ISM) in 14] as 
Q(x,t) = ( 
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where I is the 2x2 unit matrix and the 2N x 2N matrix 
S is given by 



N 



Sij = 6ijI + ^2 
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(1 < h 3 < N). Here we have introduced the following: 
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The 2x2 matrices normalized to unity in the sense of 
the square norm must take the same form as Q from their 
definition. We call them "polarization matrices" , which 
determine both the populations of the three components 
{1, 0, —1} within each soliton and the relative phases be- 
tween them. The complex constant kj denotes discrete 
eigenvalue of jth soliton, which determines a bound state 
by the potential Q in the context of ISM 0] . Cj is a real 
constant which can be used to tune the initial displace- 
ment of soliton. 

The equation @ is a completely integrable system in 
the sense that the initial value problems can be solved via 
ISM. The existence of the r matrix for this system guar- 
antees the existence of an infinite number of conservation 
laws |l it which restrict the dynamics of the system in an 
essential way. Here we show explicit forms of some con- 
served quantities: number, Nt — J dxn(x,t), n(x.i) = 
*t . $ = tr{QtQ} ; S pi nj F T = j dxf(x,t), f(x,t) = 
<j>t f<& = tv{Q* aQ}, (<x: Pauli matrices); momentum, 
P T = J dxp(x,t),p(x,t) = -ih&-d x 9 = -mtr{QtQ x } ; 
energy, E T — J dxe(x,t), e(x,t) — (h 2 /2m)d x & ■ d x <f> — 
c(n 2 + f 2 )/2 = cAy{Q\Q x - QtQQtQ}. 

If we set N = 1 in the formula (JSJ, we obtain the 
onc-soliton solution: 



Q = 2k R - 
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where e p l 2 = (2kR)~ 1 , and the subscripts R and / denote 
real and imaginary parts, respectively. We set kn > 
without loss of generality. 

In Eq. 10 , we can make out the significance of each pa- 
rameter or coordinate as follows: k^, amplitude of soli- 
ton; 2k], velocity of soliton's envelope; xr, coordinate for 
observing soliton's envelope; xi> coordinate for observing 
soliton's carrier waves. We use the term "amplitude" in 
the sense of the peak(s) height of soliton's envelope. The 
actual amplitude should be represented as ka multiplied 
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by a factor from 1 to y/2 which is determined by the type 
of polarization matrices. Note that the motion of soliton 
depends on both x and t via variables \R an d Xi: which 
elucidates the meaning of 2kj as a velocity. 

Because of the spin conservation, the one-soliton solu- 
tion can be classified by the spin states. We show that 
only two spin states are allowable, i.e., \Ft\ — Nt for 
detll = 0, and |F T | = for dctll ^ 0. 

Ferromagnetic state. — Under the condition det II = 
(a 2 = Pj), Eq. © becomes a simple form: 

Q = k R sech( X R + p/2)ne ixi - 

Now all of the mp — 0, ±1 components share the same 
wave function. Their distribution in the internal states 
reflects the elements of the polarization matrix II directly. 
One can clearly see the meaning of each parameter listed 
above. The number of particles is calculated as Nt — 
2k R . The spin of this soliton becomes 

F T = 2fc fl tr{n t crn}, |F T | = N T . (7) 

Thus, this type of soliton belongs to the ferromagnetic 
state. The momentum and the energy of the ferromag- 
netic state are = NT~hki, E r — iVrc(fcf — k a /3). 

Polar state. — In the case of det II ^ 0, a local spin 
density has one node, i.e., f(xo,t) = at a point xq — 
2M + [ln(4fcfj/| detll|) + 2e]/2k R , for each moment of t. 
Setting x' — x — xq and A^ 1 = 2\ detll|, we obtain 

f(z') = - 4fc ^ sinh(2fc * x \tr { nW|. 

[A + cosh(2fc fl a;')] 

Since each component of the local spin density is an 
odd function of x', its average value becomes zero, i.e., 
Ft = (0,0, 0) T . This implies that this type of soliton, 
on the average, belongs to the polar state jfj. Note that 
the relation Nt = 4fc# is different from that of the fer- 
romagnetic state. The momentum and the energy are 
given by P T = N T hki, E r = N T c (kj - fc|/3), respec- 
tively. The energy difference between the ferromagnetic 
state and the polar state with the same number of parti- 
cles N T is E T — E T — -N T c/l6 < 0, which is a natural 
consequence of the ferromagnetic interaction, i.e., ci < 0. 

The two-soliton solution can be obtained by setting 
N = 2 in Eq. (J5J . Since the derivation is straightforward 
but lengthy, we give an explicit formula of the general 
two-soliton solution in a separate paper 0, and here, 
focus on the two-soliton solution in the energetically fa- 
vorable ferromagnetic state (detlli = detll2 = 0), com- 
puting the asymptotic forms as t T°°i which de- 
fine a collision law of two solitons in the spinor model. 
For simplicity, we confine the spectral parameters to re- 
gions km > 0, k2R < 0, ku < 0, kii > 0, which 
correspond to a head-on collision. Under these condi- 
tions, we calculate the asymptotic forms in the final state 
(t — > oo) from those in the initial state (t — > — oo). In 



the asymptotic forms, we can consider each soliton sep- 
arately. Thus, the initial state is given by a sum of two 
solitons as Q ~ Q 1 " + Qjf, where Q 1 ? = kj R sech(xjR + 
Pj /2)H j e i ^ 1 . And for the final state, Q ~ Qf n + Qf n , 
where Qj n — fc^sech (_XjR + Pj/% + s ) H-je 1XjI . Here we 
have introduced the phase shift s = ln[l — (4fc 1 ^fc 2j R/| A;i + 
^2 1 2 ) I * r (nin| ) |] and the polarization matrices IT,- in the 
final state. Each polarization matrix IE, of the ferro- 
magnetic state can be expressed by three real variables 
Tj,0j,<Pj 13 ■ as 

n _ gir, ( cos 2 ^e-'^i cos % sin % \ 
3 \ cos ^ sin ^ sin 2 ^-e 1 ^ J ' 

In this expression, we have the following collision law: 

II, = e^Uj • uj, flj = e-'+^iij ■ uj, 

where with (j, Z) = (1,2), (2,1), 

/ cos^e" 1 ^ \ . fcj+ftf, + v 

\ sin -fe 1 2 J kj + k x 

Since each envelope is located around x ~ 2kjit, soli- 
ton 1 and soliton 2 are initially isolated at x — > ±oo and 
then travel to the opposite direction at a velocity of 2ku 
and 2fc2/, respectively. After a head-on collision, they 
pass through without changing their amplitudes and ve- 
locities and arrive at x — > ^foo in the final state. The 
collision induces rotations of their polarizations in ad- 
dition to the usual phase shifts. The collision laws for 
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FIG. 1: Time evolution of |0o| 2 (left column), \<j>i\ 2 (middle 
column), |</>-i| 2 (right column) for (a) fcj = 0.75, (b) 0.25, 
and (c) 0.05, with k R = 0.5, a 1>2 = 4/17, /3i = 72 = 16/17, 
and 71 = /? 2 = 1/17. 
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FIG. 2: to versus ki/kn for T = 1 (solid line), 0.5 (dashed 
line), 0.0157 (dash-dot line), and (dotted line). 



other cases, two-soliton of polar-ferromagnetic (detlli ^ 
0,detn 2 = 0) and polar-polar (detlli ^ 0,detn 2 ^ 0), 
can be obtained in the same way [l6| . 

Figure ^ shows the time evolution of the density pro- 
files for {1,0,-1} components in three different veloci- 
ties: (a) -k u = k 2 i = kj = 0.75, (b) kj = 0.25, and (c) 
kj = 0.05, with kiji = —k2R = ku = 0.5, ai, 2 = 4/17, 
Pi = 72 = 16/17, and 71 = /3 2 = 1/17. In the initial 
state, soliton 1 (left mover) consists mostly of the = 1 
component and, on the contrary, soliton 2 (right mover) 
almost lies mmp — ~1. A fast collision, Fig.^a), makes 
the solitons almost transparent to each other. As kj de- 
creases, the residence time inside the collisional region 
increases, and the mixing among the components occurs 
in each outgoing soliton. In Fig. ^c) , the components 
are switched between the solitons after their collision. As 
seen in Fig. \l}[h) clearly, the number of each component 
can vary not only in each soliton but also in the total 
during the collision in consequence of the spin-exchange 
interaction. This contrasts to the Manakov system [15l |. 
where the total number of each component is conserved. 

We can gain a better understanding of the two-soliton 
collision by recasting it in terms of the spin dynamics. 
The total spin conservation restricts the motion of the 
spin of each soliton on a circumference around the total 
spin axis. Since a spin of the ferromagnetic soliton is 
given by Eq. JJJ), that of the jth soliton in the initial 
state is Fj = 2 1 kj r | (sin 6j cos ipj , sin 9j sin (fj , cos 9j ) T . 
When we set \km\ = |fc 2 fl| = Nt/A, the final state spins 
Fj are obtained through Fi, 2 by Fj = cos 2 (w/2)Fj + 
sin 2 (w/2)F ; +sinw(F J x Fj)/'|F T |, where F T = Fi + F 2 , 
and w is a rotation angle of the spin precession. The ro- 
tation angle to is determined only by the ratio ki/kn and 
the magnitude of the normalized total spin T = \Ft/Nt\ 
as cos (cj/2) = e' 3 ' 2 with e s = 1 + (/c fl /A:/) 2 T 2 . 

Figure [3 shows w as a function of the ratio kj/kn for 
different values of the normalized total spin: T = 1, 0.5, 
0.0157, and 0. In consistency with Fig. ^ h exhibits 
that uj becomes larger as ki/kn decreases. The large 
(small) total spin makes the spin of each soliton rotate 
a lot (bit). When Ft is zero, corresponding to the case 



of antiparallel spin collision, the spin precession cannot 
occur as shown by the dotted line in Fig. |2] 

In this Letter, we have introduced the integrable model 
which describes the dynamics of F = 1 spinor BECs in 
one dimension. Utilizing the inverse scattering method, 
we have obtained the multiple soliton solutions. One- 
soliton solutions are classified into two distinct spin 
states: ferromagnetic, |Ft| = Nt and polar, |Ft| = 0. 
We have also shown the collision law for the two solitons 
of the ferromagnetic state and identified their collision 
with the spin precession dynamics around the total spin. 
We believe these properties should be observed in ex- 
periment and lead to a variety of applications such as 
coherent atom transport and quantum information. 

J. I. thanks T. Tsuchida for many useful discussions. 
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